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MULTIMODE  ANALYSIS  OF  BRAGG  REFLECTORS  FOR 
CYCLOTRON  MASER  APPLICATIONS 

L  INTRODUCTION 

The  cyclotron  auto-resonance  maser  (CARM)  is  a  promising  source  of  high  power 
radiation  in  the  100  GHz  to  500  GHz  frequency  range  that  may  impact  die  requirements  of 
advanced  systems  for  i^plicatitMis  such  as  radar  systems,  communications  systems  and  plasma 
heating.  The  requirements  for  guide  magnetic  field  strength  and  electron  energy  in  a  CARM  are 
advantageous  when  compared  with  competing  devices.  Compared  with  a  gyrotron,  the  required 
magnetic  field  strengdi  requirement  is  reduced  because  of  the  Dopplo*  shift  of  the  radiation. 

Because  the  CARM  depends  on  a  convective  instability,  CARM  oscillator  operation  must 
take  place  in  a  cavity  that  provides  feedback  of  the  radiation  onto  the  elecmm  beam.  For  typical 
CARM  operation,  the  restmator  reflectors  must  provide  high  reflectivity  for  modes  that  are  far 
from  cutoff.  Mode  selectivity  is  desirable,  preferably  with  disoriminatitHi  between  modes  differing 
in  either  transverse  or  axial  strocture.  Modes  that  are  near  cutoff  must  be  siq>pressed  in  order  to 
minimize  competition  from  the  gyrotron  interaction.  Finally,  the  CARM  cavity  must  allow 
unrestricted  passage  of  an  electron  beam  parallel  to  the  axis. 

CARM  resonator  reflectors  may  be  achieved  by  using  a  small  periodic  corrugation  of  the 
waveguide  surface.  Each  ripple  provides  a  small  reflection  of  wave  anqplitude.  If  the  corrugation 
period  is  such  that  the  radiation  reflected  from  each  of  the  corrugations  adds  in  phase,  the 
cmrugated  section  can  be  highly  reflective,  and  is  known  as  a  Bragg  reflector^.  This  paper  treats 
the  reflection  due  to  rippled  wall  sections  as  a  mode  conversion  from  a  forward  wave  to  a 
backward  wave,  with  an  approach  to  the  analysis  that  is  similar  to  the  analyses  of  other  mode 
conveners^.  Mode  conversion  from  the  desired  mode  to  parasitic  modes  in  the  corrugated 
section  is  included  in  the  analysis.  A  CARM  oscillator  with  a  Bragg  resonator  is  illustrated  in 
figure  1. 
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Previous  woik  on  Bragg  reflectors  has  focused  on  a  single  noiode  approach  to  the  reflection 
calculations^*  Palmer*  exanaines  the  effects  of  multiple  modes;  however,  his  derivation  requires 
fhat  the  coupled  differential  equations  for  the  mode  anqilitudes  be  solved  numerically.  This  paper 
presents  a  method  of  solving  the  coupled  noode  equations  that  only  requires  inversion  of  a  single 
4  by  4  c(Miq)lex  matrix. 


a  DERIVATION  OF  THE  COLTLED-MODE  EQUATIONS 


The  general  case  of  waveguides  with  cross  sections  that  change  along  the  axis  has  been 
treated  by  Solymar^,  who  formulated  the  equadmis  governing  the  waveguide  modes  as  a  set  of 
coupled  differential  equations  for  the  wave  amplitude,  the  “generalized  telegraphist’s" 
equations.  The  set  of  equations  consists  of  two  equations  for  each  waveguide  mode;  each 
equation  describes  the  amplitude  of  either  the  forward  or  the  backward  component  of  one 
waveguide  mode. 

where  Ai  is  the  amplimde  of  the  ith  waveguide  mode,  k  is  the  wave  number  of  the  mode,  A  is  the 
wave  impedance  of  the  mode,  and  5  is  a  wave-wave  coupling  coefficient  The  designations  and 
~  signify  forward  and  backward  going  components,  while  the  sum  p  is  over  all  waveguide  modes. 
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Equations  la  and  b  can  be  written  in  a  slowly-varying-amplitude  form  by  writing 


The  coefficient  of  the  term  due  to  the  impedance  variation  is  the  z  derivative  of  the 
logarithm  of  the  wave  inq)edance. 


^  t  tr  ^  \  t  ll\  1  ^ 


whoe  00  is  the  angular  frequency  of  the  radiation. 


If  a  sinusoidal-profile  rippled  waU  reflectOT  is  chosen. 


-a{z)~-tjkgsm{kgz-<l>) 


the  wave  impedance  term  K  becmnes 


■£ln  K  =  ^s^^sm(kgz  - 
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a-yjxi^  -  m* 


for  TE-TE  coiq)ling. 


^  -m(aifc^oSin(fcgz) 

^  acJkfk^ix^^ -m^) 


for  TE-TM  coupling,  and 


^  -kgt,(k^x/ 


for  TM-TM  coupling. 


Equaticms  2  can  be  simplified  with  the  transformation 


bi*  =  ^.c‘2* 


to  arrive  at  the  following  expression  for  the  mode  amplitudes. 


%-i^fr -o-.fr  ^iHu; 


%=-i^fr*o-,f:--2Hif;. 
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Where  G  is  the  combination  of  the  impedance  term  and  the  wall  cunwit  term,  and 

A^^2ki-k^ 


FbrTE  modes. 


G  =  y 
2 


xjl.  ~  n^a\<o^ 


h^ixl-nf) 
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and 


U  CD^Ic^+k!^ 
2a  k 
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for  the  TM  modes^.  is  the  cross  mode  coupling  tenn,  and  is 
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for  a  IE  mode  coupling  to  a  TE  mode. 


Lmkna 
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for  a  TE  mode  coupling  to  a  TM  mode,  and 
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for  a  TM  nxxie  coupling  to  a  TM  naode. 

Since  equations  10  form  a  first  order  system  of  linear  ordinary  differential  equations,  the 
solution  to  equations  10  can  be  determined  by  assuming  a  solutitm  of  the  form^ 


i«1 
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where  n  is  twice  the  number  of  modes,  y  is  an  eigenvalue  of  the  linear  system,  ^  is  the 
eigenvector  corresponding  to  y,  and  c  is  a  complex  number  needed  to  match  the  boundary 
conditions. 


In  most  cases,  only  two  modes  need  to  be  considered:  the  mode  of  interest  and  the  next 
nearest  mode  with  the  same  aamuthal  mode  number.  In  the  two  mode  case,  the  solution  of  the 
differential  equations  reduces  to  finding  the  eigenvalues  and  eigenveaors  of  the  corresponding 
matrix. 
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The  four  eigenvalues  of  this  system  are 
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y,  =  +  C?  +  G?  +  2Hi  +  -  4  -  C? + <^’)’  +  4Hi((G,  +  <^)’  -  ( A,  -  A,)’] 

Y,  =  ±i-^-A»  -  A'  +  G?  +  G|  +  2HJ.  -  ^(a’  -  ^  -  GJ  +  £?)’  +  4H|[(G,  +  G.)’  -  (A,  -  A,)'] 


.and 


and  the  corresponding  general  solution  is 


'<W' 

(a,’  -  G,"  +  y.’Xri  -  iA,) + S’(-y, + iA,)' 

«i'W 

ss  r 

[G,(A,^-G,’  +  y,’)+G,S'] 

h 

4<^(ri -<A,)+G,(y, -iA,)] 

,  S[(iy,  +  A,Xy,  +  iA,)+i(G.G,-S’)] 

C{Aj"  -  Gj"  +  Y^^Y-Y^  -  iA0+S^(y,  +  tA^) 
[Gi(A2^  +  yi^j+GjS^] 

5[G2(-ri  -  tA,)+G^(-y,  -iAj)] 

s[Hy,  +  A,X-yt  +  iA2)+i(G,G2-S^)]  J 


C3 


'  s(('■r2+^Xr,+iA,)+i{GA-s’)] 

^p2(“)'i  ~  ^i)  ■*■  ~  ^2)] 

i[(^(A,'-G,'  +  r2’)+G,S'] 

(-A,^ + G,’  -  rj'Xn + >A2)+ s'(r2 + «A,) 


\ 


s[Hy2  +  A0(-y2  +  iA2)+i(GA-S^)]  ' 

S[G2(+y2  -iA,)+G,(+y2  -iAj)] 
»|G2(a,^  - G,^  +  y2  ) + G,S^] 

(-A,"  +  G,^  -  Y2^)hy2  +  *^2) + S"(-y2  +  iA,)^ 
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The  constants  cj  through  C4  are  determined  by  the  boundary  conditions  and  are,  in 
general,  complex.  For  the  case  where  mode  1  is  incident  on  the  reflector,  the  constants  are 
determined  by  the  following  relation: 
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r«r(o)] 

0 

4(0) 

0 

<4(^)j 

where  L  is  the  length  of  the  rippled  section.  Due  to  the  complexity  of  the  solution  of  equation  19, 


constants  ci  duough  C4  are  detennined  numerically.  The  power  reflectivity  of  the  corrugated 


section  in  the  mode  of  interest  is 


gi~(0)i 


<{0) 


,  and  die  reflection  into  die  stray  mode  is 


•hm 


<(0) 


m.  EXAMPLES 


As  an  example  of  the  application  of  the  methods  in  this  paper,  we  have  calculated  the 
reflection  and  mode  conversion  for  a  TFn  mode  Bragg  reflector  at  85  GHz  and  a  TE6I  mode 
Bragg  reflector  at  100  GHz.  The  TEn  mode  is  of  interest  for  applications  where  the  election 
beam  is  centered  cm  the  axis  of  the  waveguide,  as  in  some  cyclotron  auto-resonance  masers 
(CARMs)  and  firee-electron  lasers  (FELs). 

In  the  TEii  mode  design,  the  ripple  depth  is  0.1  mm.  The  dominant  competing  mode  in  the 
reflector  section  is  the  TMn  mode,  because  it  is  the  mode  with  the  same  azimuthal  index  with 
the  closest  cutoff  frequency  to  that  of  the  TEn  mode.  The  corrugated  section  is  15  cm  long  in 
order  to  achieve  93%  reflectivity.  Table  1  contains  the  parameters  for  the  TEn  mode  reflector. 
Figure  2  illustrates  the  spatial  dependence  of  the  85  GHz  fields  in  the  reflector.  For  the  TEi  ] 
and  TMn  mode  both  the  forward  and  backward  wave  powers  are  shown.  Since  the  mode 
conversion  in  this  reflector  is  weak,  the  TMn  mode  powers  have  been  multiplied  by  ten.  The 
TEn  mode,  for  which  the  Bragg  resonance  is  satisfied,  has  a  wave  amplitude  that  decreases 
exponentially  with  distance  in  the  reflector.  The  spatial  dependance  of  the  TMn  mode  is 
oscillatory.  Figure  3  shows  the  frequency  dependence  of  the  reflection  and  transmission  of  the 
corrugated  section.  The  peak  in  the  TMn  reflected  power  occurs  where  the  Bragg  resonance  is 
satisried  for  the  TEn  to  TMn  reflection.  In  this  design  the  presence  of  the  TM  mode  has  little 


8 


effect  on  the  reflection  coefficient  Even  so,  the  conversion  of  TEi  i  to  TMi  i  provides  a  mode-mix 
of  90%  TE  mode  and  10%  TM  mode  in  the  transmitted  signal. 


Reflector  design  mode 

TEll 

Upstream  reflector  length 

4cm 

CotTugatitm  period 

1.7  mm 

Coirugatitm  anq>litude 

0.1  mm 

Center  fiequency 

85  GHz 

Power  reflectivity 

TE61 

933%% 

TM61 

03% 

Ouq)ut  mode 

TE61 

5.6% 

TM61 

0.6% 

Table  1.  Design  parameters  of  an  85  GHz  TEll  Bragg  Reflector 


The  second  example  is  a  TE6i  mode  Bragg  reflector  for  cyclotron  autoresonance  maser 
(CARM)  applications.  Cyclotron  masers  often  use  annular  electron  beams,  primarily  because 
annular  election  beams  can  be  formed  using  magnetron-injection  guns  (MlGs).  Consequently, 
resonators  designed  for  operation  in  whispering-gallery  modes  are  of  interest  for  CARMs.  Table 
2  presents  the  parameters  of  a  TE^i  mode  reflects  for  a  CARM  oscillator  resonator.  Figure  4 
shows  the  frequency  dependance  of  the  reflection  and  mode  conversion  in  the  reflector.  In  order 
to  keep  mode  conversion  to  a  tolerable  level,  the  amplitude  of  the  corrugations,  Iq  was  kept  to 
0.12  mm.  Because  the  coupling  in  the  TE^i  mode  is  strong,  the  reflector  is  comprised  of  only  24 
periods.  Consequently  the  width  of  the  resonance  is  wide  enough  to  overlap  with  the  TE^i  to 
TM^l  Bragg  resonance.  If  the  mode  conversion  is  ignored  when  the  reflectivity  is  calculated,  the 
calculated  reflectivity  is  93%.  The  resulting  mode  conversion  causes  approximately  one  half  of 
the  output  power  is  expected  to  be  in  the  TM6i  mode. 
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Figure  3  displays  a  Bragg  reflector  with  little  mode  conversion,  while  Figure  4  depicts  a 
reflector  where  the  noode  conversion  is  important  Mode  conversion  is  more  prevalent  in  the 
whispering-gallery  case  than  in  the  fundamental-mode  case  for  two  reasons.  First,  the  coupling 
coefficients  are  strongest  for  the  whispering-gallery  modes.  Second,  the  frequency  separation 
between  the  TE-TE  Bragg  resonance  and  the  nearest  TE-TM  Bragg  resonance  is  somewhat 
smaller  for  higher-order  modes.  When  the  resonance  regions  overlap,  as  they  do  in  figure  4, 
mode  conversion  can  be  large. 


Reflector  design  mode 

TEei 

Reflector  length 

4cm 

Corrugation  poiod 

1.7  mm 

Corrugation  amplitude 

0.12  mm 

Center  flequency 

100  GHz 

Power  reflectivity 

T1E61 

89.6% 

TM61 

20% 

Output  mode 

TE61 

4.0% 

TMei 

4.4% 

Table  2.  Design  parameters  of  a  TEei  Bragg  Reflector 
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IV.  CONCLUSIONS 


A  set  of  general  equations  for  the  coupled  modes  of  rippled-wall  resonators  has  been  derived. 
These  equations,  which  can  be  readily  solved  using  matrix-algebra  techniques,  demonstrate  that 
unless  they  are  designed  carefully,  Bragg  reflectors  and  resonators  can  suffer  from  mode 
conversion  that  will  reduce  the  cavity  Q-factor  for  the  desired  mode,  and  possibly  prevent 
oscillation  in  that  mode. 
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Axial  Position  (cm) 

Fig-  2  —  The  spatial  dependence  of  the  TE,,  and  TM„  electric  fields  in  an  85  GHz  Bragg  reflector  designed  to 
operate  in  the  TE„  mode.  The  tm„  mode  amplitudes  are  multiplied  by  10. 

The  dimensions  of  this  reflector  are  shown  in  Table  1 . 


/  (GHz) 

The  frequency  dependence  of  reflection  and  transmission  in  the  85  GHz  Bragg  reflector 
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